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Abstract 

Radio frequency (RF) sensing networks are a class of wireless sensor networks (WSNs) which use 
RF signals to accomplish tasks such as passive device-free localization and tracking. The algorithms 
used for these tasks usually require access to measurements of baseline received signal strength (RSS) on 
each link. However, it is often impossible to collect this calibration data (measurements collected during 
an offline calibration period when the region of interest is empty of targets). We propose adapting 
background subtraction methods from the field of computer vision to estimate baseline RSS values 
from measurements taken while the system is online and obstructions may be present. This is done by 
forming an analogy between the intensity of a background pixel in an image and the baseline RSS value 
of a WSN link and then translating the concepts of temporal similarity, spatial similarity and spatial 
ergodicity which underlie specific background subtraction algorithms to WSNs. Using experimental data, 
we show that these techniques are capable of estimating baseline RSS values with enough accuracy that 
RF tomographic tracking can be carried out in a variety of different environments without the need for 
a calibration period. 



1 Introduction 

Radio frequency (RF) tomography is a promising field of research which uses information from the strength 



of received wireless signals to construct images of physical environments 13 . RF tomography can be used 
to image static environments, but many of the technology's most interesting applications are designed to 
investigate how an environment changes over time. For instance, RF tomographic tracking can be used to 
track the movements of a person, through walls, without the need for, e.g., GPS units. 

In order for RF tomography algorithms to function, they require access to calibration values: a set of 
signal strength measurements taken when the system is deployed and transmitting, but not officially online. 
Dm'ing this time, the environment must be static and empty of temporary obstructions (i.e., in a "back- 
ground" state). Some RF tomography algorithms use these measurements directly; others compare them 
with measurements taken later, inferring how the environment is changing from how these measurements 
differ. Unfortunately, calibration values are often difficult to obtain, since it may be impossible to keep the 
environment static and empty long enough to collect them. Additionally, if the environment changes over 
time, it can be difficult to clear the area and recalibrate the system. Little research has addressed whether 
it is possible to estimate baseline RSS values without access to a calibration period. 

* A. Edelstein and M. Rabbat are with the Department of Electrical and Computer Engineering, McGill University, Montreal, 
Quebec, Canada. E-mail; andrea.edelstein@mail.mcgill.ca, michael.rabbat@mcgill.ca 
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In this paper we adapt background subtraction techniques from the field of computer vision to estimate 
RSS values on wireless links when no calibration period is available, using only measurements taken when 
the system is already online. 



1.1 Motivation 

The wireless signal measurements required for RF tomography are typically collected using wireless sensor 
networks (WSNs). These networks consist of many sensor nodes which are distributed around a region of 
interest. RF tomography applications recover important information not from any external sensor, but by 
using the nodes' own transceiver to measure the received signal strength (RSS) of the messages sent along 
wireless links between the nodes. 

Knowledge of physical laws governing how obstructions affect RSS is then exploited by RF tomographic 
imaging to build 3D representations of the surrounding area, with radio waves being used analogously to 



the X-rays employed in CT scans 19 . By comparing how RSS measurements on different links change over 



time, RF tomographic tracking follows the movements of people or objects in the region of interest 7{ [T8pO 



Although not specifically under the umbrella of RF tomography, many node localization methods also use 



RSS measurements 12 



1.2 Problem Description 

Because of the small size and relatively low cost of the nodes used in RF sensing networks, these networks 



are uniquely suited to be deployed on-the-fly 20 . For instance, they would be well-suited to supplant 
current tracking technology which may depend on GPS, RFID or video technology. While these systems 
can be accurate on the range of metres to millimetres [s], they exhibit major drawbacks: they are expensive 
and they require that tracking hardware be pre- installed. This is not the case for RF-sensing networks. 
Unfortunately, existing algorithms for RF sensing networks cannot currently work in the on-the-fly manner 
they are best-suited for because they require access to calibration data which must be collected before the 
system is brought online, during a time when no temporary obstructions are present in the network. In 
practical applications, this is often unacceptable. 

In order for RF sensing networks to achieve their true potential to be used on-the-fly, we would like to 
find a way to approximate the baseline RSS values seen when no temporary obstructions are present. 

One way to approximate this value for a given link might be to average many measurements for that link 
over time, with some no doubt taken while the RSS value is perturbed by the presence of an obstruction but 
(hopefully) many more taken while the RSS is closer to its baseline value. Over time, we might hope that 
the latter values would overwhelm the former, but this is not guaranteed. It would be better if we could 
approximate the baseline RSS for a given link by averaging only over measurements taken on that link at 
times when it is not obstructed. 

1.3 Paper Contribution and Organization 

This paper proposes and evaluates methods for online calibration of RF sensing networks. First, we consider 
a simple algorithm to estimate links' baseline RSS using the mean RSS values on each link over time. To 
improve this method's performance, we adopt background subtraction algorithms — designed to determine 
which pixels in a frame of video are currently static, background pixels — from the field of computer vision 



14 . We extend these algorithms to determine which links in a WSN are currently affected by obstructions (in 
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the "foreground" ) and which are not (in the "background" ) . By taking the mean of the measurements seen 
on a hnk only at times when it is in the background, we can achieve a better estimate of the hnk's basehne 
RSS without needing an offline cahbration period. We use experimental data to evaluate the performance 
of the proposed techniques. 

Section[2]reviews two (of many) RF sensing network applications to establish the importance of estimating 
baseline RSS. Section [3] gives an overview of background subtraction in the context of image processing, 
focusing on three particular background subtraction algorithms. Section |4] explains the modifications which 
must be made to the aforementioned algorithms to adapt them to wireless links instead of pixels. Section [s] 
presents a comparison of the performance of these background subtraction algorithms using experimental 
data. Section |6] summarizes our work and discusses future research. 



2 RF Sensing Networks 

In general, the RSS measurement vector R(fc) (in dBm) at a time k, which contains one entry for each link, 
can be decomposed into two parts, 

R^*^) = Rb - Rf^''^ (1) 

where Rb is the baseline, "background" components of the RSS signal which are static over the short term 
(although they may change over longer timescales, e.g., if the nodes shift due to wind or if their batteries 
wear down), and Rf'*^^ is the "foreground" components which are time- varying due, e.g., to people moving 



through the environment 19 . If the network is empty of temporary obstructions, Rf^*^^ equals 0. 

When a target is present, the relationship between the location of the target and its effect on Rf'-'^-' 
becomes very complicated. In general, there is no simple mathematical way to recover Rb or Rf'''^-' from 
R^'^). Most practical applications of RF sensing networks require access either to Rb or to Rf^*^-* but only 
R^'^^ can be measured directly. 



2.1 RF Tomographic Imaging and Tracking 



Wilson and Patwari's RF tomographic imaging algorithm 20 periodically inspects Rf*-'^'' to locate temporary 
obstructions such as people. This is done using an array of nodes placed around a region of interest inside 
of which objects will be moving. After the nodes are set in place, they are allowed to transmit amongst 
themselves for several minutes while no temporary obstructions are present (i.e., while Rf'*^' — 0) to estimate 
Rb. Obstructions are then allowed to enter the region of interest, while the nodes continue to probe R^'^^ 
the links' RSS values. The measured value of Rb is subtracted from R*^'^) to determine Rf^'^''. Wilson 
and Patwari assume that, for a given link (.^ any non-zero value of i?^'' [P\ is attributable to the presence of 
an obstruction located somewhere along the direct line-of-sight (LOS) path of link I. Thus, once Rf^''-' is 
known, a maximum a posteriori formulation can be used to determine the most likely obstruction location(s) 
given the observed Rf^'^'. 

Li et al.'s RF tomographic tracking method [t] operates on a similar principle. Here, successive values 
of Rf'*^' are fed to a particle filter which — given a model of the obstruction's movement — uses knowledge of 
the obstruction's prior positions to track it. 

Without a calibration period, though, and hence with no way to know the value of Rb, neither algorithm 
specifies how to isolate Rf^*^'' from R^'^). This therefore motivates our need to estimate Rb. 
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3 Background Subtraction 



Background subtraction originates from the fields of computer vision and image processing. Using video 
from a stationary camera, the goal of background subtraction is to differentiate moving "foreground" objects 
from the relatively static background. 

In its simplest form, background subtraction evaluates a test at each pixel of each frame of a video 

I^''^[n]~ B[n] (2) 

B 

where /'^'^^[ti] is the value of nth pixel in frame k (e.g., intensity on a scale from to 255 for greyscale 
video), B[n] is the known or estimated value of the true background at n and is a threshold parameter [§]. 
Evaluating this test assigns n to the foreground set T or to the background set B. 

In practice, background subtraction algorithms must account for several factors including: a background 
which is not completely static (e.g., due to fluttering leaves), noise in the image (e.g., due to camera jitter) 
and changes to the background over time (e.g., as the sun sets). Using a non-zero threshold in ([2| can 
compensate for some of these problems, but determining the correct model for B and the correct value for 
9 can be difficult. Many background subtraction algorithms have been developed to cope with these issues. 
We now describe three such methods. 

3.1 Temporal Background Modelling 

The most basic method we consider seeks only to model the distribution of the background by using a 
nonparametric Gaussian kernel density function [4] and by relying on the concept of temporal similarity, the 
assumption that if a pixel is in the background, measurements taken on it within a short time frame should 
be similar to one another. Note that this temporal similarity is not expected to hold over long periods of 
time; e.g., the background may evolve over long time-scales. More formally, this model dictates that the 
background PDF of the intensity /^'^^ [n] of a pixel n of a frame k of video can be modelled as 

PB{l^'Hn])-j^T.'^{l^'\n]-I^''%]) (3) 

i=l 

where /C(-) is a zero- mean Gaussian kernel function with a variance which is assumed to be constant. 
Here, N refers to the number of recent frames used to model the background at frame k. 
Once Pg has been formed in this way, a simple test 



PB(/W[n]) 



(4) 



is carried out to determine if each pixel in each frame is in the background. In this equation, 77 is a cost term 
which can be used to penalize or favour different classification errors, ttjt and ttb are the prior probabilities 
that a pixel is in the background or the foreground and Pjr is the foreground PDF which is assumed to be 
constant over n and k; this assumption will be modified in Sections |3.2| and |3.3| Likewise, for now, ttjt and 
TTe are considered to be constant and equal. The constants in Q can be combined into a single new constant 
9, turning the equation into 

PB{l'^''>[n])^9. (5) 
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This approach is referred to as background subtraction with temporal background modelling (TBM). 
Compared to the other background subtraction algorithms described below, TBM has lower computational 
complexity and fewer parameters which need to be tuned (TV, and 9). However, more sophisticated 
methods often yield superior performance. 



3.2 Foreground- Adaptive Background Subtraction 



Foreground-adaptive background subtraction (FABS) [8|9j uses the methods described in Section 3.1 to model 
Pg. It then models Pjr using spatial similarity, the idea that foreground pixels in the same image which 
are close together will tend to be similar in value. More specifically, McHugh et al. |9| employ the concept 
of a neighbourhood M{n): a square of a certain size centered on n, inside which all foreground pixels will 
tend to be similar. For instance, say we know that pixel n is part of the face of a person in the foreground. 
Any foreground pixels in the neighbourhood of n will then likely be flesh-coloured as well. Thus, spatial 
similarity is the assumption that foreground objects occupy space (a group of adjacent pixels) and are 
spatially homogenous in value. 

To model the foreground, the TBM test ([s]) is first applied to a given frame to assign a preliminary label 
to each pixel. Then, an initial model of the foreground, Pj^^ , is built according to 



P^. (l^"^ W) "° (6) 



where A/jto {n) is a set containing all the neighbours of n which have been assigned a foreground label by 
the preliminary TBM test. Again, /C(-) is a zero-mean Gaussian kernel function, this time with constant 
variance cr^. Once Pjr^ is calculated, it is used in 

P^,{m[n])^'' '^'> 

to assign new labels to each pixel. These labels are then used to build Afj^^, Pj^n A/jpg, PjTa, ... until a 
convergence criterion is satisfied. McHugh [sj suggests that only 2-3 iterations of this process are usually 
needed for the labels to stabilize to what can be called Pjr (with no iterative subscript). 

While the sets of background and foreground pixels returned by FABS should be closer to the ground truth 
than those returned by TBM, note that FABS introduces positive feedback to the background subtraction 
process. FABS is specifically designed to look for foreground pixels which have erroneously been left in the 
background by a simpler background subtraction algorithm. These pixels are then added to J- . This means 
that if 7] is set too high, FABS may eventually place every pixel in the image in J- . The threshold 77 must be 
chosen carefully to avoid this. 

We point out here that FABS is more computationally-complex than TBM. It also introduces several 
more parameters which must be set, namely 77, and neighbourhood size/structure. 



3.3 Foreground- Adaptive Background Subtraction with Markov Modelling of 
Change Labels 

Foreground-adaptive background subtraction with Markov modelling of change labels (FABS-MMCL) [8j|9] 
employs the aforementioned methods to model Pg and Pjr and then models ttjt and tt^ by relying on spatial 
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ergodicity, the idea that pixels which are close together will tend to belong either to the foreground or the 
background together. 

McHugh et al. [oj propose that the FABS test in Equation ([t]) first be used to label each pixel in a given 
frame. McHugh [s] then points out that, by assuming the pixels' true labels obey the Markov property, these 
labels represent a realization of a Markov random field G. The a priori probabilities of such a field taking a 
certain value follow the Gibbs distribution 

P(G = 5) = ^exp|— (8) 

where Z is a normalization constant and 7 is the distribution's natural temperature. Q corresponds to a 
pixel neighbourhood, and V (•) is an indicator function which equals if pixels n and n' have the same label 
and 1 if they do not. This prior penalizes pixels which have labels that are very different from their neighbors 
(i.e., showing a high degree of non-ergodicity) , and discourages isolated false positives. 

By temporarily assuming that G is known for all pixels except n, ([s]) can be combined with Q to write 



ryexp 



(^(K^W|-K^W|))- (9) 



Here \Mj^{n)\ and |A/B(n)| are the number of foreground and background neighbours of n as labelled by the 
first apphcation of the FABS test in Q. Labels are then updated via the FABS-MMCL test 

McHugh et al. |9j recommend repeatedly applying the test (|9| to refine the labels (a technique known as 
iterated conditional modes [2]). Although the recommended number of iterations is relatively low (« 10), 
FABS-MMCL is the most computationally-intensive of the three methods and has the most parameters to 
set. 

Note that FABS-MMCL is designed in part to prune back the effects of the overly-greedy FABS algorithm. 
Because of this, the values of 77, al and the neighbourhood size which produce the best results for FABS in 
isolation may not be optimal when FABS outputs are used as input to FABS-MMCL. 



4 Adaptation of Background Subtraction to WSNs 

Before adapting background subtraction to WSNs, we ask whether any simpler techniques can be used to 
determine Rb when it cannot be measured directly. 

One naive approach would be simply to average over all available RSS measurements for a given link. 
However, if an obstruction is moving around or across the link, this average will incorporate the obstruction's 
effects. We will still introduce this basic mean approximation (MA) approach as a starting point since it 
is easy to implement and nonetheless an improvement over the current situation wherein no ways have yet 
been proposed to determine Rb at all. 

Another naive approach would examine the RSS for a given link over a window of time, searching for 
periods when the RSS drops suddenly. One might then attribute these drops to the attenuation caused by 
a newly-appearing obstruction, use a threshold to discard these values and take the mean of the remaining 
measurements. Unfortunately, the effect of obstructions on RSS is not that straightforward; due to multi- 
path effects which can cause constructive interference, obstructions can cause RSS to decrease or increase. 
Decreasing battery voltage and shifts in the nodes' positions will also cause changes to the RSS, but we do 
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not wish to discard measurements which are perturbed for these reasons. Finally, a considerable amount 
of noise is seen even on unobstructed links [s]. These factors drastically reduce the performance of simple 
techniques such as thresholding 10 . 



4.1 Basic Background Subtraction in WSNs 

The naive algorithms mentioned above are not always equipped to deal with the complex behaviour of 
wireless signals. However, we can see strong parallels between the problem of determining which links in a 
network are currently affected by a moving obstruction (i.e., are in the foreground) and that of determining 
which pixels in an image are currently displaying a moving foreground object, pointing to the utility of 
background subtraction in the WSN domain. Before background subtraction can be applied to a WSN, it 
is necessary to establish what constitutes a "pixel," a measurement of "intensity" and a "frame" in this 
context. We can make the association between a single pixel n in an image and a single bidirectional link 
£ in a WSN, and we can see that the discretized value of the RSS (in dBm) measured on a link during an 
as-yet-to-be-defined frame k, R'-''^[£], represents an analogue to the greyscale intensity of a pixel, /'^'^-'[n]. 

In the image processing domains, one frame of video comprises one measurement for each pixel in the 
video image. Ideally, we would like our corresponding RSS frame to contain one RSS measurement for each 
bidirectional link with all these measurements taken at exactly the same time, giving us an instantaneous 
capture of the RSS values on all the links in the network. Because WSNs are subject to interference and 
packet collisions when multiple nodes transmit at the same time, this is not realistic. 

In practice, RSS measurements are usually obtained as follows: for a network with M nodes, node 1 
broadcasts a packet to all other nodes. This establishes measurements for the RSS on M— 1 unidirectional 
links stretching from node 1 to node 2, from node 1 to node 3, etc. Then, node 2 broadcasts a packet 
to all the other nodes, yielding another series of M— 1 unidirectional RSS measurements. This continues 
until all M nodes have broadcast, at which point node 1 broadcasts again. Therefore, it is impossible to 
take simultaneous measurements on all the links in the network. Fortunately, WSNs are capable of taking 
many such measurements per second: in our experiments, a 22-node network could complete such a cycle 
in « 120 ms. Because this time span is so short, we assume the conditions in the network (i.e., obstruction 
locations) do not change in the time it takes to complete one measurement cycle. 

We must also address that, because WSNs are subject to dropped packets, a single measurement cycle 
may not return a value for every link. If this happens, we do not want to wait for additional cycles to 
complete our frame. Thus, instead of a frame being comprised of one measurement for each link, we define 
a frame as being built from M consecutive measurement vectors, where M is the number of nodes in the 
network. Ideally, each node contributes one measurement vector to each frame. However, if the message from 
one node is dropped then that node may be omitted and another node included multiple times. In this case, 
we average the multiple measurement vectors for the node which appears in the same window twice. Then, 
the two directed measurements for each transceiver pair are averaged together to create one measurement for 
each uni-directional link. Finally, if measurements have been dropped for one or more links, these values are 
copied from the previous frame. The next frame is then created from the next M consecutive measurement 
vectors and so on. A consequence of this method is that there is no overlap in the measurements used from 
frame to frame, except where dropped packets force us to reuse values. 

Having established these analogies for pixels, pixel intensities and frames, we see that the temporal 
similarity underlying TBM should still hold here to some degree: if the effects of noise are small, the 
consecutive measurements taken on a given link should not vary while that link stays in the background 
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(discounting long-term changes such as those caused by a dying battery). In that case, we can calculate the 
background PDF Pg via an updated version of (|3| such that 

i=l 

and we can build our sets of background and foreground links via an updated version of ([s]) such that 

PB(i?W[^]) (11) 

Recall that N is the number of recent frames used to model the background at frame k and R^''''[i] is the 
measured RSS value for link £ at frame k. 

4.2 Spatial Similarity in WSNs 

While determining analogues for pixels, pixel intensity and frames is relatively straight-forward, applying 
spatial similarity and the concept of neighbourhoods to WSNs is more complicated: there is no reason why 
two foreground links which are close together must have similar RSS values. Consider a network of nodes 
deployed around the perimeter of a square. Let £i be the link connecting the nodes at (0,0) and (0, 1), let 
£2 be the link connecting the nodes at (0, 0) and (0, 7) and let both these links be obstructed by the same 
object. The two overlapping links are physically close together, but because is 1 m long while £2 is 7 m 
long, the log-distance path loss model suggests that their RSS values will be quite different. 

Instead, to capture the notion of similarity between links, we define the neighbourhood £{£) of link £ to 
be the set of links similar in length to £. Letting represent the length of link £ (i.e., distance between 
transmitter and receiver), C{£) is thus defined as 

C{£)^{£' ■.\di~di,\<T}, (12) 

where r > is a parameter to be tuned. This definition implies the creation of £^-'(£), the set of foreground 
neighbours associated with C{£) in frame k. In this way, spatial similarity in the image domain becomes 
similarity-of-length in the WSN domain. 

When formulating £{£) in this way, we are aware that not all links in this type of neighbourhood are 
guaranteed to have similar RSS values. Even for links of equal length, if one 1 m link in the network passes 
through a wall while another passes through free space, their values will be different. Keep in mind, though, 
that even in the image processing domain, not all foreground pixels in a neighbourhood are guaranteed to 
have a similar intensity (e.g., a flesh-coloured pixel may be next to a red pixel from a shirt). 

Furthermore, we are aware that two links of the same length whose LOS paths are both obstructed by 
the same moving object may have different RSS values if the object intersects them in different ways, e.g.. 



if an obstruction is located at the end of one link and in the center of the other 15 . In practice, though, we 
find that this difference is within the range of the normal variations due to noise seen on the links. 

We can now present new expressions for the model of the preliminary foreground PDF Pjr^ (originally in 
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Figure 1: A graphical representation of FABS-MMCL-R. The set of rectangles p associated with each 
link is shaded. We see \pi\ = 2, \p2\ = 7 and |pi 0/921 = 1, creating shared rectangle percentage 
Xieij2) = = 0.5. Similarly, jpsl = 8, |p2 Hpal = 4 and x(^2,4) = t^tJ^ = 0.57. 



min(7,8) 



6])) and for the FABS test (originally in ^) such that 



P 



and 



(13) 



(14) 



4.3 Spatial Ergodicity in WSNs 

We also consider the concept of spatial ergodicity which underlies FABS-MMCL. Spatial ergodicity carries 
over to WSNs more directly than does spatial similarity since we still say that if two links are "close together" , 
they are likely to pass through or avoid the same obstructions and hence have correlated foreground or 
background states. Nonetheless, we need to reconsider what it means for two links to be "close together." 
Must they share a common endpoint? Must their LOS paths overlap or intersect? The following subsections 
describe two methods for determining whether two links are close enough together for spatial ergodicity to 
hold. 



4.3.1 Rectangle-Based Neighbourhoods 

The first method involves partitioning the interior of the region of interest into rectangles whose side lengths 
are determined by the node separations. These rectangles are similar to the cubic voxels proposed by Wilson 



and Patwari in 20 . We denote by p the set of rectangles through which a link I passes. Then, for two 



links, ii and £2 with associated sets of rectangles pi and p2, we define x (^ij ^2), the percentage of rectangles 
common to both links as 

|pi np2| 



x(^i,4) 



X 100. 



(15) 



mindpil , \p2\) 

If this percentage is larger than a certain tuneable threshold, we deem €1 and £2 to be neighbours. A graphical 
representation of this process can be seen in Fig. [Ij We refer to this method as FABS-MMCL-R. 

Unlike in the computer vision and image processing domains, the WSN-domain neighbourhoods used for 



Figure 2: A possible relationship between an obstruction and two different links. Links £i and £2 may pass 
through a large enough percentage of common rectangles to qualify as neighbours. Nonetheless, as seen at 
the top of the figure, the links' paths may eventually diverge enough that obstructions (like the one pictured) 
do not affect them both. 



FABS and FABS-MMCL are not the same sets. As previously mentioned, we denote the neighbourhood for 
FABS as £(£). We now denote the neighbourhood used for FABS-MMCL as S{e). For FABS-MMCL-R, this 
neighbourhood is defined as 

S{1)^{£' ■.x{i,n>C} (16) 

where < C < 100 is left as a parameter to be tuned. We also define S^\£) and sj^\£), the sets of 
foreground and background neighbours, respectively, associated with S{£) in frame k. This leads to the 
FABS-MMCL test ^ becoming 



J" 

^ r/exp 



(l(|4^)(£)|-|5f(^)|)) 



(17) 



4.3.2 Obstruction-Based Neighbourhoods 



This section describes an alternative, data-driven way to define the neighborhood set S{£) used by FABS- 
MMCL. Although FABS-MMCL-R is fairly intuitive, it may have problems. For instance, if two LOS links 
have a small angle of separation, they may be close enough to cross the same rectangles over a large fraction 



of their length, thereby qualifying as neighbours according to ( 16 ). However, this is misleading if their paths 



diverge significantly by the time one of them crosses the spot where an obstruction is located, as seen in 
Fig. [2] This section discusses an alternative approach to determine neighborhoods based on estimates of 
where the obstruction may be located. 

Recall that we do not have any a priori information about the location of these obstructions. We can still 
hypothesize, though, that their locations may be inferred by looking for rectangles which contain particularly 



large numbers of foreground links. Specifically, just before each application of ( 17), we can count the number 



of foreground links passing through each rectangle in p. We can then find the rectangle r in p with the largest 
number of foreground links passing through it and assign S^^\£) to be the set of links which pass through 
r. This is shown graphically in Fig. [3] We refer to this method as FABS-MMCL-0. 

Of course, r may not be located near an obstruction at all since, in fact, it is quite possible that none 
of the rectangles in p are near the obstruction (e.g., if the LOS of £ is positioned entirely on the extreme 
east side of the network while the obstruction is on the west side of the network). If this is the case, even 



though r may contain more foreground links than any other rectangle in p. 



is still likely to be small. 

This will cause the exponent term in the FABS-MMCL test (seen in (17)) to be highly negative, decreasing 



SP{£) 



the right hand side of the equation. In this way, no harm is done if r is artificially deemed to contain an 
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Figure 3: A graphical representation of FABS-MMCL-0. Link £ is the vertical link. Other foreground links 
(as designated by FABS or by previous iterations of FABS-MMCL-0) are also shown. Excluding £ itself, 
one pixel in p is crossed by 3 foreground links, one pixel in p is crossed by 1 foreground link and the rest 
of the pixels in p are crossed by foreground links. Since max(3, 1, 0) = 3, r is the rectangle crossed by 3 
foreground links. 



obstruction. However, in our experiments with the TBM method, we notice that even when we know for 
a fact that r contains an obstruction, and even when our classification of foreground and background links 
results in an accurate estimate of the baseline RSS, there may still be relatively few foreground links and 
hence many background links passing through r. This means that, if we apply FABS-MMCL-0 to the output 



of the TBM algorithm as described, the exponential term in (17) would almost always be highly- negative 



To counteract this tendency, we introduce a term to the FABS-MMCL test such that 



^ Tyexp 
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SP{£) 



(18) 



Recall that B and F are the sets of background and foreground links respectively. Then, by resetting 
p—\B\ — \F\ before each application of (18 1, we ensure that the distribution of the exponent is roughly 



centered around instead of always being negative. 



4.3.3 Other Methods of Building Neighbourhoods 

Although we present two models for building S{1), we are aware that other researchers have considered 
related questions. In [6], Gudmundson proposes a model for quantifying the correlation of the shadowing 
observed on different links, but this model only applies to links with at least one common endpoint, which 
is not the case for most links in a WSN. Wang et al. 17 extend Gudmundson's work to model shadowing 
correlation between two links, £i and ^2, when £2 is created by moving the nodes which form £1 to a new 
location. We might apply this model to a case where £1 and £2 are two links defined by four completely 
separate senders and receivers, but Wang et al.'s model also assumes that the distance between £1 and £2 is 
small compared with the lengths of the links themselves. Again, this does not apply to most links in a WSN. 
Agrawal and Patwari (l 11 also propose the network shadowing (NeSh) model to calculate the shadowing 
correlation seen across links in a wireless network, and this may also be used as an alternative method to 
define link neighborhoods. 
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5 Experimental Evaluation 



5.1 Experiment Setup 

Through an association with the Beijing University of Posts and Telecommunications (BUPT) in Beijing, 
China, we have access to data from a number of RF sensing experiment sets carried out by BUPT students. 
Each set of experiments used the same methodology, hardware and software. Only the number of nodes and 
the environment in which they were deployed varied between experiment sets. 

The specific hardware employed in these experiments consisted of a fieet of Texas Instruments CC2530 
System-on-Chips equipped with built-in RSS indicators capable of measuring RSS to a resolution of 1 dBm. 
These nodes conform to the IEEE 802.15.4 standard, transmitting in the 2.4 GHz band with a configurable 



output power 16 . For these experiments, the nodes were set to transmit at dBm and each node was 
augmented with an external antenna. 

The nodes were programmed to follow the token-ring protocol from Section [4.1[ and all M nodes commu- 
nicated with one another, forming a fully-connected topology which could directly measure the RSS between 
every pair of nodes in the network. The system was configured to probe this RSS (by having the next node 
in the ring send a new wireless message) approximately every 5 ms. The nodes also sent their measurements 
to an M-|-lst node, which forwarded the data via USB cable to a laptop, which recorded the values. This 
extra node acted only as an intermediary to the laptop; any RSS values on links involving this node are 
irrelevant to the experiments. 

These experiments were conducted as follows: In each new environment, the nodes were mounted 
on w 1 m-high stands around the perimeter of the region of interest. They were then allowed to trans- 
mit amongst themselves for several minutes while the network was kept vacant, measuring the baseline RSS 
on each link directly. After this offline calibration period was complete, a person entered the region of inter- 
est and walked around the nodes. This person followed a path marked on the ground, and a timestamp was 
recorded as they passed various points on their route so that their location was known at all times. Each 
experiment set thus consists of baseline RSS data collected while the region of interest was empty, coupled 
with data collected during one or more "walks." Note that the length of each walk varied; some lasted only 
long enough for « 150 frames of measurements to be obtained, while in other cases, up to 700 frames of 
measurements were collected. 

The first set of experiments was conducted in an outdoor field at BUPT. Outdoor environments such as 
fields are often the most conducive to RF tomography applications as the lack of walls and furniture means 
that troublesome multipath effects are far less prominent. In this case, 22 nodes were placed in a 7 m x 7 m 
square, as seen in Fig. |6] The odd node positioning in the upper-right corner is due to the fact that two 
nodes were discovered to be defective after the experiments had been run and the measurements taken by 
these nodes had to be discarded. The walker's path was a 6 m x 6 m square inside the nodes. 

The second set of experiments was conducted in an indoor lab, with 24 nodes positioned in a 3 m x 3.6 m 
rectangle (see Fig.|9]). The walker's path was a 1.8 m x 2.4 m rectangle inside the nodes. Although there was 
furniture along the walls of the room (outside the region being monitored by the nodes) , no furniture was 
located inside the region of interest. Thus, the only thing obstructing the LOS paths of the wireless signals 
was the walker sent into the region of interest; nonetheless — as in indoor environments in general — the walls, 
floor, ceiling and furniture still contribute to the multipath effects which generally make RF sensing more 
challenging indoors. 

The third set of experiments was also conducted in an indoor lab. This time, some of the nodes were 
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Figure 4: Node deployment for experiment set 3. 



placed in a hallway outside the lab (see Fig. |4| , thus monitoring the inside of the lab through a solid wall 
which passed through the gap between nodes 6 and 7 and 1 and 22 in Fig. |4(c) This scenario is one of the 
hardest for RF sensing applications to cope with. 

The data from these experiments allow us to test algorithms in several different environments. 



5.2 Performance Evaluation Methods 

One way to evaluate our algorithms would be to compare the links they assign to the foreground set J- with 
ground truth knowledge of which links are actually in the foreground. However, there is no way to know this. 
One might assume that links whose LOS paths intersect the obstruction ought to be in the foreground, but 
this is not guaranteed, especially indoors where multipath effects are non-trivial and the non-line-of-sight 
(NLOS) components of a signal can be heavy contributors to its RSS value. If one of these NLOS components 
is obstructed, the associated link may well be in the foreground. We therefore evaluate our algorithms in 
two ways. We first look at their ability to estimate the baseline RSS measured during the offline calibration 
period. Then, to see how useful our values are to a real-world application, we examinine the performance of 
an RF tomographic tracking algorithm initialized with our values. 

5.2.1 Metric 1: Comparison With Measured Rb 

Because we have access to direct measurements of the baseline RSS from the calibration period of each 
experiment set, we define -R_b[^] as the mean RSS measured on £ over this time period, when the region of 
interest was intentionally left vacant. Recall that this is the data which existing RF tomographic algorithms 
already assume is available to them. 

Our estimate of the baseline RSS on link I, i?_B[£], is then calculated according to 

where K is the total number of frames available at a given time and V'*-'''' [^] is an indicator function 

,1, if£e^ in frame 
V'('=)[£] = <| , (20) 

0, if ^ e in frame k 
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built using the output of our background subtraction algorithms. This ip'-'^'^ll] will take different values 
depending on whether background subtraction is carried out using TBM, FABS or FABS-MMCL. 

Once we obtain Rb[^], we can calculate the root mean square estimation error of the approximation 
(measured in dBm/link) such that 



Estimation Error ^ J - RbI^])^ (21) 



V eeA 

where A is the total number of links in the network. This provides us with a numerical qualification for 
how well we have estimated the baseline RSS, which can easily be compared across our different algorithms. 
We can also calculate this metric for the mean approximation (MA) method described in Section [4j which 
constitutes the simplest possible way to estimate baseline RSS without access to a calibration period. Note 
that the MA algorithm corresponds to setting ip^'^^ [i] = 1 for all values of k and i. 

5.2.2 Metric 2: Tracking Performance 

To get an intuitive feel for the implications of a certain level of estimation error, we use our data in a real- 
world RF sensing network application. Many such applications exist, but we choose to employ Li et al.'s 
RF tomographic tracking algorithm 7 , which must be given a value for the baseline RSS on each link in 
the network before it can begin tracking. By comparing tracking performance when using Rb, the vector 
of measured baseline RSS values, to tracking performance when using different versions of Rb, we can see 
whether our estimates of the baseline RSS are accurate enough for real-world applications. 

Because Li et al.'s algorithm is non-deterministic, for each value of Rb or Rb, we run A = 100 realizations 
of the algorithm and report the performance over this ensemble in the form of the root mean square tracking 
error calculated as 

A 



Tracking Error 



a=l \ t=l 



1^ 



(22) 



where zj"'' is the estimated value of the person's position at time f e T as reported by realization a. 
5.2.3 Additional Evaluation 

To further understand the behaviour of our background subtraction algorithms, we look at which links and 
RSS values are selected for the foreground using histograms and other diagrams. We also investigate the 



algorithms' run-times in Section 5.4 



Because this work constitutes the first examination of the potential of background subtraction to estimate 
baseline RSS in WSNs, we are mostly interested in determining the maximum improvements which can 
be gained from this technique. Therefore, we present results which are obtained when the background 
subtraction parameters {6, N, etc.) are specifically those found (via a series of parameter sweeps) to give 
the best performance and lowest estimation error. To offset this, in Section [5. 5[ we provide an examination 
of the algorithms' sensitivity to parameter values which differ from these settings. 
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Figure 5: RMSE obtained when estimating the basehne RSS for experiment set 1. The 143 frames used were 
collected while a walker made one round trip around the set path over a period of ~ 20 seconds. 



Table 1: Parameter values used for experiment set 1. 



Parameter 


TBM 


FABS 


FABS- 
MMCL-R 


N 


25 


25 


25 




4 


4 


4 


9 


0.17 


0.17 


0.17 


T 














4 


4 


1] 




1 


1 


7 






5 


C 






95 



5.3 Results for Each Experiment Sets 
5.3.1 Results for Experiment Set 1 

Fig. [5] presents the estimation error obtained when using our background subtraction algorithms to estimate 
baseline RSS values for an experiment conducted in an outdoor field. Table [T] presents the parameter values 
used to obtain the pictured results. 

From Fig. [5] we see that MA is fairly accurate already and that TBM provides another sizable improve- 
ment compared to MA. Recall that these are the least computationally-complex of all the algorithms we 
present and that no knowledge of the nodes' positions is required to apply them. Consequently, it is partic- 
ularly encouraging that they are this successful. We also see that the more complicated algorithms do not 
actually do any better than TBM in this case. 

Since this original experiment set includes only a single "walk" , we validate these results by repeating 
our analysis on data sets collected several months later, in a different outdoor field also located at BUPT; 
these results are shown in j^. Note that in both fields, we see the same pattern in our results: TBM shows 
an improvement over MA while the more complex FABS-based approaches do not show any significant 
improvement compared to TBM. Additionally, it is worth noting that these similar results are obtained in 
two different outdoor environments with nodes set out in two slightly different patterns (since the broken 
nodes 11 and 13 from experiment set 1 were fixed before more experiments were conducted), but using the 
same parameters seen in Table [T] This is very encouraging, because even if these parameter values seem 
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Table 2: RMSE for tracking using experiment set 1. 



Algorithm 


Tracking Error (m) 


Calibration Data 


0.4073±0.0655 


MA 


0.4363±0.0711 


TBM 


0.4319±0.0655 


FABS 


0.4335±0.0667 


FABS-MMCL-R 


0.4339±0.0616 




02468 02468 02468 

(m) (m) (m) 



(a) TBM foreground (b) FABS fore- (c) FABS-MMCL-R 
ground foreground 




02468 02468 02468 

(m) (m) (m) 



(d) TBM back- (e) FABS back- (f) FABS-MMCL-R 
ground ground background 

Figure 6: Links assigned to the foreground and background by various background subtraction algorithms 
for a single representative frame in experiment set 1. The known position of the obstruction during this 
frame is marked. 

slightly arbitrary, this shows that at least they are transferrable between similar environments. 

Returning to the experiment set 1 data, we also use Li et al.'s tracking algorithm [T] (initialized in turn 
with Rb and with various values of Rb) to track a walking person. Tracking errors are shown in Table [2j 

From these results, we can see that the ideal situation (where we initialize tracking with measured 
data) provides the best performance, as expected. We can also see that using MA to build Rb causes our 
performance to suffer only slightly, while TBM provides a small improvement over MA — though even this 
small improvement may well be desirable in systems intended to replace video or RFID tracking systems. 
Finally, we see that the more complex algorithms are no more helpful than TBM. In all cases, the tracking 
error is low, lending support to the idea that any of the four methods shown can serve as an excellent 
substitute to measuring Rb directly. This makes sense in an outdoor environment; the situation will change 
in more complex, indoor environments. 

To understand why FABS and FABS-MMCL do not improve over MA and TBM, we examine which 
links are being selected as foreground and background links by each algorithm. This is seen (for a single 
representative frame) in Fig. |6] Note that any link which is not in the foreground in a given frame must be 
in the background, so the bottom row of figures is the logical negation of the top row. 

As we described in Section |3.3[ FABS is designed to grow the set of foreground links created by TBM. 
However, in Fig.|6j we see that, in this case, FABS is removing links from the set created by TBM. We might 
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Figure 7: Histogram of RSS measurements seen on the link connecting nodes at (7,1) and (1,0) during 
experiment set 1. 
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Figure 8: RMSE obtained when estimating the basehne RSS for experiment set 2. The 693 frames used were 
collected while a walker made four round trips around the set path over a period of ~ 80 seconds. 



hypothesize that this is occurring because TBM is being overly-greedy, perhaps due to 6 being set too high. 
Recall though, that the parameters have all been adjusted to optimize performance. Looking at Fig. |6]just 
provides more supporting evidence, then, that for outdoor environments, the low-complexity TBM algorithm 
is already doing an admirable job identifying foreground links, and any attempt to improve on its efforts 
using FABS or FABS-MMCL is futile. 

Having examined the background vs. foreground behaviour of all the links in the network for a single 
frame, we will also look at the behaviour of one link in the network over all frames. In Fig. 7(a) we see 



a histogram of the RSS measured on a particular link in the WSN (in this case, the link connecting nodes 
at (7, 1) and (1, 0)) over all frames in and in Fig. |7(b)[ we see a histogram of the subset of these values 
deemed to be in the background by TBM. From these figures, we see that the outliers from Fig. 7(a)| are 
properly identified as being in the foreground and are removed in Fig. |7(b)| This behaviour is representative 
of that seen on all links in the network. 



5.3.2 Results for Experiment Set 2 

Fig. [8] presents the root mean square estimation error obtained when using our background subtraction 
algorithms to estimate baseline RSS values for an experiment conducted in an indoor lab. Table [3] presents 
the associated parameter values. 

If we compare Fig. [8] to the similar figure generated for experiment set 1, we can see some important 
differences. First, if we compare the performance of MA in each case, we see that even this simple approach 
performs much better in the outdoor environment. This hints that our task of properly estimating the 
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Table 3: Parameter values used for experiment set 2. 



Parameter 


TBM 


FABS 


FABS- 


FABS- 








MMCL-R 


MMCL-0 


N 


35 


35 


35 


35 




17 


17 


17 


17 





0.05 


0.05 


0.05 


0.05 


T 




0.75 


0.75 


0.75 






10 


10 


10 


V 




5 


5 


1 


7 






25 


50 


C 






35 





baseline RSS is far harder indoors than it is outdoors. Fortunately, TBM provides an improvement over MA 
in the indoor case and, moreover, FABS and FABS-MMCL provide further improvement still (with FABS- 
MMCL-0 providing more improvement than FABS-MMCL-R). Even after applying these algorithms, the 
lowest estimation error we achieve is still higher than the estimation error seen for the outdoor experiment 
sets, but this is understandable given the increased complexity and more intense multipath effects inherent 
in an indoor environment; additional results showing similar behavior are presented in 

We can see another result of these multipath effects if we compare the convergence rates between the 
indoor and outdoor experiment sets. This convergence rate depends on the amount of noise affecting the 
RSS measurements and, in particular, on how often we have access to background measurements for each 
link. At one extreme, if there is absolutely no measurement noise and if no links are ever deemed to be in the 
foreground, Rb can converge after a single frame. On the other hand, if many links are in the foreground in 
each frame (e.g., because there are many obstructions, because the obstruction(s) are very large or because 
multipath effects cause the obstruct ion (s) to affect many links) or if the same links are in the foreground 
for many frames (again, because the obstruction(s) affect many links or because they do not move around 
very much), it may take several frames before we can obtain a background measurement for each link. The 
estimate Rb will then converge more slowly. 

In our case, the same size obstruction was moving at approximately the same rate in both the indoor and 
outdoor sets of experiments. However, due to multipath effects, more links are in the foreground for each 
frame in the indoor environment (as will be seen when we graphically examine the foreground links selected 
for a single frame). Therefore, fewer measurements are available for use in calculating Rb at any given time, 
causing slower convergence. 

Again, we also use Li et al.'s tracking algorithm to detect the movements of a person walking around 
the lab. The tracking errors seen over 100 realizations of the algorithm (initialized with Rb and Rb) are 
shown in Tabled 

As expected, using MA in place of the measured Rb increases the tracking error from 0.6060 m to 
0.6801 m. Simple TBM reduces this error to 0.6404 m while the more complicated FABS-based approaches 
all reduce the error closer to the level of (or below the level of!) the error seen when Rb is used. The error 
variance is also reduced using more complex methods. 

Once again, we also examine which links are assigned to the foreground and background by each algorithm. 
This is presented in Fig. |9] 
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Table 4: RMSE for tracking using experiment set 2. 



Algorithm 


Tracking Error (m) 


Calibration Data 


0.6060±0.0022 


MA 


0.6801±0.1763 


TBM 


0.6404±0.1962 


FABS 


0.5994±0.0033 


FABS-MMCL-R 


0.6221±0.0841 


FABS-MMCL-0 


0.6120±0.0956 




024 024 024 

(m) (m) (m) 



(a) TBM foreground (b) FABS fore- (c) FABS-MMCL-O 
ground foreground 




024 024 02 

(m) (m) (m) 



(d) TBM back- (e) FABS back- (f) FABS-MMCL-O 
ground ground background 

Figure 9: Links assigned to the foreground and background by various background subtraction algorithms 
for a single representative frame in experiment set 2. The known position of the obstruction during this 
frame is marked. 
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Figure 10: RMSE obtained when estimating the basehne RSS for experiment set 3. The 388 frames used 
were collected while a walker made two round trips around the path in Fig. 4(c) over a period of ~ 50 
seconds. 



Table 5: Parameter values used for experiment set 3. 



Parameter 


TBM 


FABS 


FABS- 


FABS- 








MMCL-R 


MMCL-0 


N 


35 


35 


35 


35 




1 


1 


1 


1 





0.05 


0.05 


0.05 


0.05 


T 




1 


1 









1 


5 


10 


V 




5 


15 


4 


7 






50 


10 


C 






15 





Unlike with our outdoor experiment set, in this case, we can see that the background subtraction algo- 
rithms are behaving as intended. Namely, TBM selects an initial set of foreground links, FABS expands this 
set, and FABS-MMCL then prunes the set back. Although it may seem like we end up with an unintuitively- 



large proportion of links being assigned to the foreground in Fig. 9(c) we can't argue with the fact that 
large foreground sets such as this provide the closest approximation to the measured baseline RSS. Note, 
however, that this successful approximation is not just a matter of carelessly assigning a large proportion of 
the network's links to the foreground. For instance, no matter how we change the parameters of TBM to 
increase the number of links it assigned to the foreground, we could not decrease the algorithm's estimation 
error to the level achieved by FABS-MMCL. The finesse of the more complex approaches is required in order 
to improve the baseline RSS estimate in this indoor environment. 



5.3.3 Results for Experiment Set 3 

Fig.[lO]presents the estimation error obtained when using our background subtraction algorithms to estimate 
the baseline RSS values for experiments in which measurements were collected indoors, through a solid wall. 
Table [5] presents the associated parameter values. 

In this case, we have once again increased the complexity of our environment, which also increases the 
lowest achievable estimation error seen in Fig. [T0| Nonetheless, basic TBM still provides an improvement 
over MA. Likewise, as in the non-through-wall indoor experiment, we can obtain further improvements over 
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Table 6: Lost track percentage and RMSE for tracking using experiment set 3. 



Algorithm 


Lost Track % 


Tracking Error (m) 


Calibration Data 


36 


2.7606±0.0116 


MA 


100 


N/A 


TBM 


95 


1.8364±0.8932 


FABS 


93 


2.6290±0.0081 


FABS-MMCL-R 


74 


1.3068±0.2080 


FABS-MMCL-0 


29 


1.7380±0.8578 




(111) (m) (m) 

(a) TBM foreground (b) FABS fore- (c) FABS-MMCL-O 
ground foreground 




(111) (111) (111) 



(d) TBM back- (e) FABS back- (f) FABS-MMCL-O 
ground ground background 



Figure 11: Links selected as being in the foreground and background by various background subtraction 
algorithms for a single representative frame in experiment set 3. The position of the obstruction during this 
frame is marked. 



TBM by using FABS and FABS-MMCL-O. 

Tracking also becomes more complicated with through- wall measurements. It now takes the particle filter 
a long time to converge to the proper track, if indeed it ever does. For the data we collected, whether we 
initialize tracking with Rb ot with Rb itself, the algorithm inevitably requires ~ 300 timesteps before the 
particle filter finds the correct track. Furthermore, over our 100 realizations initialized with Rb, in 36 cases, 
the algorithm never converged at all, resulting in very obviously "lost" tracks. 

In this environment, it thus becomes important to consider not just the tracking error for different Rb's, 
but also the percentage of lost tracks. This is summarized in Tabic [6] For instance, when Rb is calculated 
using MA, every single one of the 100 realizations is a lost track. Similarly, when Rb is calculated using 
TBM or with FABS, 95% and 93% of the realizations are lost tracks; consequently, the associated tracking 
error values are not reliable. It is only the values of Rb produced by FABS-MMCL that allow the tracking 
algorithm to converge in >25% of realizations. Then, in these cases, the tracking error is actually consistently 
lower than the error seen when the tracking algorithm is initialized using Rb. 

In Fig. 11 we present a one- frame snapshot of which links in the through- wall environment are being 
assigned to the foreground and background by each algorithm. Comparing this to the similar figure presented 
for the indoor environment, we can see that the layout of the foreground is similar in both cases. 



21 



Table 7: Observed per-frame run times for each algorithm 



Exp. 


TBM 


FABS 


FABS-MMCL- 


FABS-MMCL- 


Set 


(ms) 


(ms) 


R (ms) 


(ms) 


Set 1 


1.056 


43.3561 


240.1544 




Set 2 


1.3952 


23.7115 


160.0984 


432.9784 


Set 3 


1.278 


62.7504 


415.7762 


441.8479 



5.4 Algorithm Speed 

As previously mentioned, our various background subtraction algorithms differ substantially in terms of the 
amount of computation they require. Table [7] presents the average runtime for each algorithm for a single 
frame (measured over 50 repetitions run on a dual-core 2.13 GHz Intel processor with 2 GB of RAM). An 
analysis as to how these runtimes depend on the algorithm parameters can be found in [sj. 

5.5 Sensitivity to Parameter Values 

The results reported above show the performance of the proposed background subtraction methods when 
parameter settings are tuned to yield the best performance. Additional experiments considering the effects 
of parameter mismatch and investigating sensitivity to parameter values are reported |3]; they are omitted 
here due to space limitations. The main conclusion of these experiments is that the proposed methods are 
most sensitive to changes in the decision thresholds {6, rj), the natural temperature 7 used in the Gibbs 
distributions, and the size of neighborhoods used in FABS-MMCL-R. 

6 Conclusion 

In this work, we explained how several pre-existing background subtraction algorithms could be adapted for 
use in determining the background, baseline RSS values of links in a WSN while the network is already online 
and obstructions may already be present. This constitutes the first set of methods proposed to estimate 
these values when they cannot be measured during an offline calibration period. Using experimental data, 
we have demonstrated that background subtraction techniques can successfully estimate the baseline RSS 
in a range of different environments. 

We can also make recommendations as to which of our proposed algorithms is best-suited to a given 
environment and application. Since MA and TBM are both at least moderately effective at estimating the 
baseline RSS in experiment sets 1 and 2, and since these are the least computationally-complex algorithms — 
in addition to being the only ones suited for applications such as node localization — these are probably the 
best choice for non-through- wall applications; if more accurate RSS measurements are needed, e.g., for use 
with a non-robust RF-sensing application, TBM, FABS or FABS-MMCL may be preferred. For through- 
wall applications, the complexity of FABS-MMCL — and in particular of FABS-MMCL-0 — will be required. 
Thus, all of the methods we proposed fill their own particular niches. 

In addition to examining the performance of our methods, we also obtained a rough estimate of their 
convergence periods. This will vary depending on what is occurring inside the region of interest, but in all 
our experiments, it took less than 300 frames (« 40 seconds) for our estimates to converge. Waiting 40 
seconds to obtain estimates for the baseline RSS represents a huge improvement over having to wait for the 
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entire network to be evacuated to obtain these values (if evacuating the network is even possible at all). 

Finally, we presented some representative runtimes for our different algorithms, and — from our pre- 
liminary experiments with different parameter values — we started to gain some insight into how sensitive 
background subtraction for WSNs is with respect to each algorithm parameter. 



6.1 Future Work 

In the future, we would like to explore more complex ways of modelling the spatial similarity and spatial 
ergodicity between links, in the hope of seeing better performance in indoor environments. For instance, we 
are currently experimenting with creating weighted neighbourhoods which preserve the idea that some of a 
link's neighbours may be more similar to it (in length, in physical location, in covariance as predicted by the 



NeSh model 11 , etc.) than others. 

There is also work to do to further our understanding of how to optimize the background subtraction 
parameters. Although we conducted a brief study of these parameters, a true understanding of their effects 
will not be complete without an experimental survey which looks at the ideal parameters for different types 
of environments; for different types, sizes and numbers of obstructions; for different network layouts, etc. 

Even in the absence of these improvements and extensions, however, we can still say that using background 
subtraction to estimate the baseline RSS in wireless sensor networks which are already online is an interesting 
new technique which has already yielded promising results and which fills a definite need in the domain of 
RF sensing network applications. 
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